A Chromosome-Level Genome Assembly of the Non-Hematophagous Leech Whitmania pigra (Whitman 1884): Identification and Expression Analysis of Antithrombotic Genes.

Despite being a non-hematophagous leech, Whitmania pigra is widely used in traditional Chinese medicine for the treatment of antithrombotic diseases. In this study, we provide a high quality genome of W. pigra and based on which, we performed a systematic identification of the potential antithrombotic genes and their corresponding proteins. We identified twenty antithrombotic gene families including thirteen coagulation inhibitors, three platelet aggregation inhibitors, three fibrinolysis enhancers, and one tissue penetration enhancer. Unexpectedly, a total of 79 antithrombotic genes were identified, more than a typical blood-feeding Hirudinaria manillensis, which had only 72 antithrombotic genes. In addition, combining with the RNA-seq data of W. pigra and H. manillensis, we calculated the expression levels of antithrombotic genes of the two species. Five and four gene families had significantly higher and lower expression levels in W. pigra than in H. manillensis, respectively. These results showed that the number and expression level of antithrombotic genes of a non-hematophagous leech are not always less than those of a hematophagous leech. Our study provides the most comprehensive collection of antithrombotic biomacromolecules from a non-hematophagous leech to date and will significantly enhance the investigation and utilization of leech derivatives in thrombosis therapy research and pharmaceutical applications.

Observation of High-Capacity Monoclinic B-Nb2 O5 with Ultrafast Lithium Storage.

Apart from Li4 Ti5 O12 , there are few anode substitutes that can be used in commercial high-power lithium-ion batteries. Orthorhombic T-Nb2 O5 has recently been proven to be another substitute anode. However, monoclinic B-Nb2 O5 of same chemistry is essentially inert for lithium storage, but the underlying reasons are unclear. In order to activate the "inert" B-Nb2 O5 , herein, nanoporous pseudocrystals to achieve a larger specific capacity of 243 mAh g-1 than Li4 Ti5 O12 (theoretical capacity: 175 mAh g-1 ) are proposed. These pseudocrystals are rationally synthesized via a "shape-keep" topological microcorrosion process from LiNbO3 precursor. Compared to pristine B-Nb2 O5 , experimental investigations reveal that B-Nb2 O5- x delivers ≈3000 times higher electronic conductivity and tenfold enhanced Li+ diffusion coefficient. An ≈30% reduction of energy barrier for Li-ion migration is also confirmed by the theoretical calculations. The nanoporous B-Nb2 O5- x delivers unique ion/electron transport channels to proliferate the reversible and deeper lithiation, which activate the "inert" B-Nb2 O5 . The capacitive-like behavior is observed to endow B-Nb2 O5- x ultrafast lithium storage ability, harvesting 136 mAh g-1 at 100 C and 72 mAh g-1 even at 250 C, superior to Li4 Ti5 O12 . Pouch-type full cells exhibit the energy density of ≈251 Wh kg-1 and ultrahigh power density up to ≈35 kW kg-1 .

Mycolicibacterium arseniciresistens sp. nov., isolated from lead-zinc mine tailing, and reclassification of two Mycobacterium species as Mycolicibacterium palauense comb. nov. and Mycolicibacterium grossiae comb. nov.

A Gram-positive, acid-fast, aerobic, rapidly growing and non-motile strain was isolated from lead-zinc mine tailing sampled in Lanping, Yunnan province, Southwest China. 16S rRNA gene sequence analysis showed that the most closely related species of strain KC 300T was Mycolicibacterium litorale CGMCC 4.5724T (98.47 %). Additionally, phylogenomic and specific conserved signature indel analysis revealed that strain KC 300T should be a member of genus Mycolicibacterium, and Mycobacterium palauense CECT 8779T and Mycobacterium grossiae DSM 104744T should also members of genus Mycolicibacterium. The genome size of strain KC 300T was 6.2 Mb with an in silico DNA G+C content of 69.2 mol%. Chemotaxonomic characteristics of strain KC 300T were also consistent with the genus Mycolicibacterium. The average nucleotide identity, digital DNA-DNA hybridization and average amino acid identity values, as well as phenotypic, physiological and biochemical characteristics, support that strain KC 300T represents a new species within the genus Mycolicibacterium, for which the name Mycolicibacterium arseniciresistens sp. nov. is proposed, with the type strain KC 300T (=CGMCC 1.19494T=JCM 35915T). In addition, we reclassified Mycobacterium palauense and Mycobacterium grossiae as Mycolicibacterium palauense comb. nov. and Mycolicibacterium grossiae comb. nov., respectively.

The Coulomb effect in nonsequential double ionization by counter-rotating two-color elliptical polarization fields.

Using a three-dimensional classical ensemble model, nonsequential double ionization (NSDI) of Ar atoms by counter-rotating two-color elliptical polarization (TCEP) fields is investigated. The major axes of the two elliptical fields are aligned in different directions. The relative alignment of the two elliptical fields strongly affects the waveform of the combined electric field and the ultrafast dynamics of NSDI in TCEP fields. Numerical results show that the correlated electron momentum distributions in the x direction evolve from a V-shaped structure near the axis to a distribution concentrated on the diagonal with the angle between the two elliptical major axes increasing. The asymmetry of the energy sharing between the two electrons during recollision results in the V-shaped structure in the correlated momentum spectrum. Back analysis indicates that the recollision times of a part of the trajectories move from the peak to the valley of the combined electric field with the angle between the two elliptical major axes increasing. Therefore, for the case of a larger angle between the two elliptical major axes, the electrons experience a longer time to escape away from the vicinity of the parent ion and thus the stronger Coulomb effect from the parent ion makes the momentum difference between two electrons small, which results in a distribution concentrated on the diagonal. This provides an effective avenue to control the electron ultrafast dynamics in NSDI.

Revisiting the Asian Buffalo Leech (Hirudinaria manillensis) Genome: Focus on Antithrombotic Genes and Their Corresponding Proteins.

Leeches are well-known annelids due to their obligate blood-feeding habits. Some leech species secrete various biologically active substances which have important medical and pharmaceutical value in antithrombotic treatments. In this study, we provided a high-quality genome of the Asian buffalo leech (Hirudinaria manillensis), based on which we performed a systematic identification of potential antithrombotic genes and their corresponding proteins. Combining automatic and manual prediction, we identified 21 antithrombotic gene families including fourteen coagulation inhibitors, three platelet aggregation inhibitors, three fibrinolysis enhancers, and one tissue penetration enhancer. A total of 72 antithrombotic genes, including two pseudogenes, were identified, including most of their corresponding proteins forming three or more disulfide bonds. Three protein families (LDTI, antistasin, and granulin) had internal tandem repeats containing 6, 10, and 12 conserved cysteines, respectively. We also measured the anticoagulant activities of the five identified hirudins (hirudin_Hman1 ~ hirudin_Hman5). The results showed that three (hirudin_Hman1, hirudin_Hman2, and hirudin_Hman5), but not the remaining two, exhibited anticoagulant activities. Our study provides the most comprehensive collection of antithrombotic biomacromolecules from a leech to date. These results will greatly facilitate the research and application of leech derivatives for medical and pharmaceutical purposes in the treatment of thrombotic diseases.

Chryseobacterium luquanense sp. nov., a casein-hydrolysing bacterium from the Jiaozi Mountain in Yunnan, PR China.

Strain KC 927T was isolated during an investigation of the soil bacteria diversity on Jiaozi Mountain, central Yunnan, Southwest China. The strain was Gram-stain-negative, rod-shaped, non-motile, oxidase-negative, catalase-positive and aerobic. Results of 16S rRNA gene alignment and phylogenetic analysis indicated that strain KC 927T was a member of the genus Chryseobacterium and closely related to Chryseobacterium caseinilyticum GCR10T (98.4%), Chryseobacterium piscicola DSM 21068T (98.3 %) and 'Chryseobacterium formosus' CCTCC AB 2015118T (97.9 %). With a genome size of 4 348 708 bp, strain KC 927T had 33.5 mol% DNA G+C content and contained 4012 protein-coding genes and 77 RNA genes. The average nucleotide identity and digital DNA-DNA hybridization values between strain KC 927T and C. caseinilyticum GCR10T, C. piscicola DSM 21068T and 'C. formosus' CCTCC AB 2015118T were 80.1, 79.6 and 90.7 %, and 25.5, 23.6 and 42.0 %, respectively. The main polar lipid of strain KC 927T was phosphatidylethanolamine and the respiratory quinone was MK-6. The major fatty acids (≥10 %) were iso-C15 : 0, iso-C17 : 1 ω9c and iso-C17 : 0 3-OH. Evidence from phenotypic, phylogenetic and chemotaxonomic analyses support that strain KC 927T represents a new species of the genus Chryseobacterium, for which the name Chryseobacterium luquanense sp. nov. is proposed. The type strain is KC 927T (=CGMCC 1.18760T=JCM 35707T).

Chromosome-level Dinobdella ferox genome provided a molecular model for its specific parasitism.

BACKGROUND: Dinobdella ferox is the most frequently reported leech species parasitizing the mammalian nasal cavity. However, the molecular mechanism of this special parasitic behavior has remained largely unknown. METHODS: PacBio long-read sequencing, next-generation sequencing (NGS), and Hi-C sequencing were employed in this study to generate a novel genome of D. ferox, which was annotated with strong certainty using bioinformatics methods. The phylogenetic and genomic alterations of D. ferox were then studied extensively alongside the genomes of other closely related species. The obligatory parasitism mechanism of D. ferox was investigated using RNA-seq and proteomics data. RESULTS: PacBio long-read sequencing and NGS yielded an assembly of 228 Mb and contig N50 of 2.16 Mb. Along Hi-C sequencing, 96% of the sequences were anchored to nine linkage groups and a high-quality chromosome-level genome was generated. The completed genome included 19,242 protein-coding genes. For elucidating the molecular mechanism of nasal parasitism, transcriptome data were acquired from the digestive tract and front/rear ends of D. ferox. Examining secretory proteins in D. ferox saliva helped to identify intimate connections between these proteins and membrane proteins in nasal epithelial cells. These interacting proteins played important roles in extracellular matrix (ECM)-receptor interaction, tight junction, focal adhesion, and adherens junction. The interaction between D. ferox and mammalian nasal epithelial cells included three major steps of pattern recognition, mucin connection and breakdown, and repair of ECM. The remodeling of ECM between epithelial cells of the nasal mucosa and epithelial cells of D. ferox may produce a stable adhesion environment for parasitism. CONCLUSIONS: Our study represents the first-ever attempt to propose a molecular model for specific parasitism. This molecular model may serve as a practical reference for parasitism models of other species and a theoretical foundation for a molecular process of parasitism.

Computational immune synapse analysis reveals T-cell interactions in distinct tumor microenvironments.

The tumor microenvironment (TME) and the cellular interactions within it can be critical to tumor progression and treatment response. Although technologies to generate multiplex images of the TME are advancing, the many ways in which TME imaging data can be mined to elucidate cellular interactions are only beginning to be realized. Here, we present a novel approach for multipronged computational immune synapse analysis (CISA) that reveals T-cell synaptic interactions from multiplex images. CISA enables automated discovery and quantification of immune synapse interactions based on the localization of proteins on cell membranes. We first demonstrate the ability of CISA to detect T-cell:APC (antigen presenting cell) synaptic interactions in two independent human melanoma imaging mass cytometry (IMC) tissue microarray datasets. We then generate melanoma histocytometry whole slide images and verify that CISA can detect similar interactions across data modalities. Interestingly, CISA histoctyometry analysis also reveals that T-cell:macrophage synapse formation is associated with T-cell proliferation. We next show the generality of CISA by extending it to breast cancer IMC images, finding that CISA quantifications of T-cell:B-cell synapses are predictive of improved patient survival. Our work demonstrates the biological and clinical significance of spatially resolving cell-cell synaptic interactions in the TME and provides a robust method to do so across imaging modalities and cancer types.

Preparation of fish decalcified bone matrix and its bone repair effect in rats.

Decalcified bone matrix has great potential and application prospects in the repair of bone defects due to its good biocompatibility and osteogenic activity. In order to verify whether fish decalcified bone matrix (FDBM) has similar structure and efficacy, this study used the principle of HCl decalcification to prepare the FDBM by using fresh halibut bone as the raw material, and then degreasing, decalcifying, dehydrating and freeze-drying it. Its physicochemical properties were analyzed by scanning electron microscopy and other methods, and then its biocompatibility was tested by in vitro and in vivo experiments. At the same time, an animal model of femoral defect in rats was established, and commercially available bovine decalcified bone matrix (BDBM) was used as the control group, and the area of femoral defect in rats was filled with the two materials respectively. The changes in the implant material and the repair of the defect area were observed by various aspects such as imaging and histology, and its osteoinductive repair capacity and degradation properties were studied. The experiments showed that the FDBM is a form of biomaterial with high bone repair capacity and lower economic cost than other related materials such as bovine decalcified bone matrix. FDBM is simpler to extract and the raw materials are more abundant, which can greatly improve the utilization of marine resources. Our results show that FDBM not only has a good repair effect on bone defects, but also has good physicochemical properties, biosafety and cell adhesion, and is a promising medical biomaterial for the treatment of bone defects, which can basically meet the clinical requirements for bone tissue repair engineering materials.

Nuclear localization of mitochondrial TCA cycle enzymes modulates pluripotency via histone acetylation.

Pluripotent stem cells hold great promise in regenerative medicine and developmental biology studies. Mitochondrial metabolites, including tricarboxylic acid (TCA) cycle intermediates, have been reported to play critical roles in pluripotency. Here we show that TCA cycle enzymes including Pdha1, Pcb, Aco2, Cs, Idh3a, Ogdh, Sdha and Mdh2 are translocated to the nucleus during somatic cell reprogramming, primed-to-naive transition and totipotency acquisition. The nuclear-localized TCA cycle enzymes Pdha1, Pcb, Aco2, Cs, Idh3a promote somatic cell reprogramming and primed-to-naive transition. In addition, nuclear-localized TCA cycle enzymes, particularly nuclear-targeted Pdha1, facilitate the 2-cell program in pluripotent stem cells. Mechanistically, nuclear Pdha1 increases the acetyl-CoA and metabolite pool in the nucleus, leading to chromatin remodeling at pluripotency genes by enhancing histone H3 acetylation. Our results reveal an important role of mitochondrial TCA cycle enzymes in the epigenetic regulation of pluripotency that constitutes a mitochondria-to-nucleus retrograde signaling mode in different states of pluripotent acquisition.

Analysis on the Vulnerability of a Tunnel Entrance under Internal Explosion.

Tunnels play an essential role in the transportation network. Tunnel entrances are usually buried at a shallow depth. In the event of an internal explosion, the blast pressure will cause severe damage or even collapse of the tunnel entrance, paralyzing the traffic system. Therefore, an accurate assessment of the damage level of tunnel entrances under internal blast loading can provide effective assistance for the anti-blast design of tunnels, post-disaster emergency response, and economic damage assessment. In this paper, four tunnel entrance specimens were designed and fabricated with a scale ratio of 1/5.5, and a series of field blast tests were carried out to examine the damage pattern of the tunnel entrances under internal explosion. Subsequently, static loading tests were conducted to obtain the maximum bearing capacity of the intact specimen and residual bearing capacities of the post-blast specimens. After that, an explicit non-linear analysis was carried out and a numerical finite element (FE) model of the tunnel entrance under internal blast loading was established by adopting the arbitrary Lagrangian-Eulerian (ALE) method and validated based on the data obtained from the field blast and static loading tests. A probabilistic vulnerability analysis of a typical tunnel entrance subjected to stochastic internal explosions (assuming various charge weights and detonation points) was then carried out with the validated FE model. For the purpose of damage assessment, the residual bearing capacity of the tunnel entrance was taken as the damage criterion. The vulnerability curves corresponding to various damage levels were further developed based on the stochastic data from the probabilistic vulnerability analysis. When the charge weight was 200 kg, the tunnel entrance exhibited slight or moderate damage, while the tunnel entrance suffered severe or even complete damage as the charge weight increased to 1000 kg. However, the tunnel entrance's probability of complete damage was less than 10% when the TNT charge weight did not exceed 1000 kg.

Electron correlation and recollision dynamics in nonsequential double ionization by counter-rotating two-color elliptically polarized laser fields.

Nonsequential double ionization (NSDI) of Argon atoms by counter-rotating two-color elliptically polarized (TCEP) fields is investigated with a three-dimensional classical ensemble model. Different from two-color circularly polarized fields, the combined electric field in TCEP pulses has no symmetry and the ionized electron mainly returns to the parent ion from one direction. Thus the electron momentum distributions show strong asymmetry. Numerical results show with the increase of the relative phase between the two elliptical fields, the return angle of the travelling electron, i.e., the angle between the return direction of the electron and the +x direction, gradually decreases. Moreover, the dominant behavior of electron pairs evolves from anti-correlation to correlation with the relative phase increasing. This provides an avenue to control the return angle and electron correlation behavior by the relative phase between the two elliptical fields.

Siccirubricoccus soli sp. nov., a novel bacterial species isolated from Jiaozi Mountain soil.

An isolate of Gram-stain-negative and strictly aerobic bacterium, designated KC 17139T, was isolated from Jiaozi Mountain sample in Yunnan, China. Cells were non-motile cocci to oval, catalase-positive and oxidase-positive. Growth occurred at 0-7% NaCl (w/v; optimum, 0%), pH 6.0-8.0 (optimum, pH 7.0) and 15-45 °C (optimum, 28-37 °C). The polar lipids were diphosphatidylglycerol (DPG), phosphatidylethanolamine (PE), phosphatidylglycerol (PG), phosphatidylcholine (PC) and four unidentified aminolipids (UALs). Strain KC 17139T contained summed feature 8 (comprising C18:1 ω7c and/or C18:1 ω6c), C18:1 2OH and C16:0 as major cellular fatty acids (> 5%) and ubiquinone-10 as the sole isoprenoid quinone. The 16S rRNA gene sequence analysis indicated that strain KC 17139T shared highest similarities with Siccirubricoccus phaeus 1-3T (96.7%) and Siccirubricoccus deserti SYSU D8009T (95.0%). Strain KC 17139T clustered with the two Siccirubricoccus type strains, but formed a separate branch in both 16S rRNA gene and genome-scale phylogenetic dendrograms. The genomic DNA G + C content of strain KC 17139T was 71.2%. Genomic comparisons between strain KC 17139T and its close relatives showed the highest digital DNA-DNA hybridisation to S. phaeus (35.5%), highest average nucleotide identity to S. phaeus (88.2%), indicating that strain KC 17139T represents a novel species. On the basis of results of phenotypic, chemotaxonomic and molecular analysis, we report a new bacterium strain KC 17139T belonged to genus Siccirubricoccus, for which the name Siccirubricoccus soli sp. nov. is proposed. The type strain is KC 17139T (= CGMCC 1.18756T = JCM 35132T).

A new species of medicinal leech in the genus Hirudo Linnaeus, 1758 (Hirudiniformes, Hirudinidae) from Tianjin City, China.

Medicinal leeches in the genus Hirudo have been utilized for therapeutic procedures for thousands of years. To date, six known species of Hirudo are widely distributed in different regions of the Eurasian continent. In this study, a new medicinal leech species Hirudotianjinensis Liu, sp. nov. is described based upon specimens collected from Tianjin City, China. The new species can be distinguished from its congeners by a combination of characters: blackish green dorsum with five continuous yellow longitudinal stripes; six sensillae on dorsal annulus a2 of segments VIII-XXV; greyish green ventrum with irregular bilateral dark brown spots; dorsum and abdomen separated by a pair of pale yellow stripes; front half atrium wrapped by white prostate; apparent albumen gland; epididymis massive in relation to ejaculatory bulb. The phylogenetic tree based upon COI implies a sister relationship to H.nipponia Whitman, 1886. A key to the known species is provided.

Plin2-mediated lipid droplet mobilization accelerates exit from pluripotency by lipidomic remodeling and histone acetylation.

Metabolic switch is critical for cell fate determination through metabolic functions, epigenetic modifications, and gene expression. However, the mechanisms underlying these alterations and their functional roles remain unclear. Here, we show that Plin2-mediated moderate lipid hydrolysis is critical for pluripotency of embryonic stem cells (ESCs). Upon exit from pluripotency, lipid droplet (LD)-associated protein Plin2 is recognized by Hsc70 and degraded via chaperone-mediated autophagy to facilitate LD mobilization. Enhancing lipid hydrolysis by Plin2 knockout promotes pluripotency exit, which is recovered by ATGL inhibition. Mechanistically, excessive lipid hydrolysis induces a dramatic lipidomic remodeling characterized by decreased cardiolipin and phosphatidylethanolamine, which triggers defects in mitochondrial cristae and fatty acid oxidation, resulting in reduced acetyl-CoA and histone acetylation. Our results reveal how LD mobilization is regulated and its critical role in ESC pluripotency, and indicate the mechanism linking LD homeostasis to mitochondrial remodeling and epigenetic regulation, which might shed light on development and diseases.

Insights into gut microbiota communities of Poecilobdella manillensis, a prevalent Asian medicinal leech.

AIMS: Medicinal leeches (Annelida: Hirudinea) are fresh water ectoparasitic species which have been applied as traditional therapy. However, gut microbiota could bring high risks of opportunistic infections after leeching and arouses great interests. Here, gut bacterial and fungal communities of an Asian prevalent leech Poecilobdella manillensis were characterized and analysed through culture-independent sequencing. METHODS AND RESULTS: With high coverage in 18 samples (>0.999), a more complicated community was apparent after comparing with previous leech studies. A total of 779/939 OTUs of bacteria and fungi were detected from leech guts. The bacterial community was dominated by the phylum Bacteroidetes and Synergistetes. Genera Mucinivorans and Fretibacterium accounted mostly at the genus level, and genus Aeromonas showed an extremely low abundance (2.02%) on average. The fungal community was dominated by the phylum Ascomycota and Basidiomycota. At the genus level, the dominant OTUs included Mortierella, Geminibasidium and Fusarium. The analysis of core taxa included those above dominant genera and some low-abundance genera (>1%). The functional annotation of the bacterial community showed a close correlation with metabolism (34.8 ± 0.6%). Some fungal species were predicted as opportunistic human pathogens including Fusarium and Chaetomiaceae. CONCLUSIONS: The present study provides fundamental rationales for further studies of such issues as bacteria-fungi-host interactions, host fitness, potential pathogens, and infecting risks after leeching. It shall facilitate in-depth explorations on the safe utilization of leech therapy. SIGNIFICANCE AND IMPACT OF STUDY: Present paper is the first-ever exploration on microbiota of a prevalent Asian medicinal leech based on culture-independent technical. And it is also the first report of gut fungi community of medicinal leech. The diversity and composition of bacteria in P. manillensis was far different from that of the European leech. The main components and core OTUs indicate a particular gut environment of medicinal leech. Unknown bacterial and fungal species were also recovered from leech gut.

The genome of medicinal leech (Whitmania pigra) and comparative genomic study for exploration of bioactive ingredients.

BACKGROUND: Leeches are classic annelids that have a huge diversity and are closely related to people, especially medicinal leeches. Medicinal leeches have been widely utilized in medicine based on the pharmacological activities of their bioactive ingredients. Comparative genomic study of these leeches enables us to understand the difference among medicinal leeches and other leeches and facilitates the discovery of bioactive ingredients. RESULTS: In this study, we reported the genome of Whitmania pigra and compared it with Hirudo medicinalis and Helobdella robusta. The assembled genome size of W. pigra is 177 Mbp, close to the estimated genome size. Approximately about 23% of the genome was repetitive. A total of 26,743 protein-coding genes were subsequently predicted. W. pigra have 12346 (46%) and 10295 (38%) orthologous genes with H. medicinalis and H. robusta, respectively. About 20 and 24% genes in W. pigra showed syntenic arrangement with H. medicinalis and H. robusta, respectively, revealed by gene synteny analysis. Furthermore, W. pigra, H. medicinalis and H. robusta expanded different gene families enriched in different biological processes. By inspecting genome distribution and gene structure of hirudin, we identified a new hirudin gene g17108 (hirudin_2) with different cysteine patterns. Finally, we systematically explored and compared the active substances in the genomes of three leech species. The results showed that W. pigra and H. medicinalis exceed H. robusta in both kinds and gene number of active molecules. CONCLUSIONS: This study reported the genome of W. pigra and compared it with other two leeches, which provides an important genome resource and new insight into the exploration and development of bioactive molecules of medicinal leeches.

MAP2K6 remodels chromatin and facilitates reprogramming by activating Gatad2b-phosphorylation dependent heterochromatin loosening.

Somatic cell reprogramming is an ideal model for studying epigenetic regulation as it undergoes dramatic chromatin remodeling. However, a role for phosphorylation signaling in chromatin protein modifications for reprogramming remains unclear. Here, we identified mitogen-activated protein kinase kinase 6 (Mkk6) as a chromatin relaxer and found that it could significantly enhance reprogramming. The function of Mkk6 in heterochromatin loosening and reprogramming requires its kinase activity but does not depend on its best-known target, P38. We identified Gatad2b as a novel target of Mkk6 phosphorylation that acts downstream to elevate histone acetylation levels and loosen heterochromatin. As a result, Mkk6 over-expression facilitates binding of Sox2 and Klf4 to their targets and promotes pluripotency gene expression during reprogramming. Our studies not only reveal an Mkk phosphorylation mediated modulation of chromatin status in reprogramming, but also provide new rationales to further investigate and improve the cell fate determination processes.